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Symmetry protected topological (SPT) phases are gapped quantum phases with a certain symme- 
try, which can all be smoothly connected to the same trivial product state if we break the symmetry. 
For non-interacting fermion systems with time reversal (T), charge conjugation (C), and/or U(l) 
(N) symmetries, the total symmetry group can depend on the relations between those symmetry 
operations, such as TNT' 1 = N or TNT' 1 = — N. As a result, the SPT phases of those fermion 
systems with different symmetry groups have different classifications. In this paper, we use Kitaev's 
K-theory approach to classify the gapped free fermion phases for those possible symmetry groups. 
In particular, we can view the U(l) as a spin rotation. We find that superconductors with the 
Sz spin rotation symmetry are classified by Z in even dimensions, while superconductors with the 
time reversal plus the S z spin rotation symmetries are classified by Z in odd dimensions. We show 
that all 10 classes of gapped free fermion phases can be realized by electron systems with certain 
symmetries. We also point out that to properly describe the symmetry of a fermionic system, we 
need to specify its full symmetry group that includes the fermion number parity transformation 
(— ) JV . The full symmetry group is actually a projective symmetry group. 

PACS numbers: 71.23.-k, 02.40.Re 
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I. INTRODUCTION 

We used to believe that all possible phases and phase 
transitions are described by Landau symmetry breaking 
theory. 1-3 However, the experimental discovery of frac- 
tional quantum Hall states 4 ' 5 and the theoretical discov- 
ery of chiral spin liquids 6, indicate that new states of 
quantum matter without symmetry breaking and with- 
out long range order can exist. Such new kind of orders is 
called topological order, 8,9 since their low energy effective 
theories are topological quantum field theories. 10 At first, 
the theory of topological order was developed based on 



their robust ground state degeneracy on compact spaces 
and the associated robust non-Abelian Berry's phases. 8,9 
Later, it was realized that topological order can be char- 
acterized by the boundary excitations, 11,12 which can be 
directly probed by experiments. One can develop a the- 
ory of topological order based on the boundary theory. 13 

Since its introduction, we have been trying to ob- 
tain a systematic understanding of topological orders. 
Some progresses have been made for certain simple cases. 
We found that all 2D Abelian topological orders can 
be classified by integer if-matrices. 14-16 The 2D non- 
chiral topological orders (which can be smoothly con- 
nected to time reversal and parity symmetric states) are 
classified by spherical fusion category. 17-20 The recent 
realization of the relation between topological order and 
long range entanglement 19 (defined through local uni- 
tary transformations 21,22 ) allows us to separate another 
simple class of gapped quantum phases - symmetry pro- 
tected topological (SPT) phases. SPT phases are gapped 
quantum phases with a certain symmetry, which can 
all be smoothly connected to the same trivial product 
state if we break the symmetry. A generic construc- 
tion of bosonic SPT phases in any dimension using the 
group cohomology of the symmetry group was obtained 
in Ref. 23,24. The constructed SPT phases include in- 
teracting bosonic topological insulators and topological 
superconductors (and much more). 

Another type of simple systems are free fermion sys- 
tems, for which a classification of gapped quantum 
phases can be obtained through K-theory 25-27 or non- 
linear er-model of disordered fermions. 28 They include 
the non-interacting topological insulators 29-35 and the 
non-interacting topological superconductors. 36-40 Most 
gapped quantum phases of free fermion systems are SPT 
phases protected by some symmetries, such as topolog- 
ical insulators protected by the time reversal symme- 
try. While some others have intrinsic topological orders 
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Insulators with spin-orbital coupling and spin order (or non-coplanar spin order) 
(lc-m • (tc 4 / + icjn 2 ■ o-Cj' + icj.n 3 • ctcv) + (c-nj • ad + cjn 2 ■ crcj + c' k nz ■ erc k ) 


17(1) 


Superconductors with spin-orbital coupling and spin order (or non-coplanar spin order) 
c «i • o-Ci + c n 2 • crcj + c! n 3 • <xc fc + (c-|-iC W - cjiCfj J 


"none" = Z^ 


T T T 

Insulators with spin-orbital coupling iqrii • (rev + ic'n2 • fc,/ + icj.n.3 • erc ft / 
(symmetry: charge conservation and time reversal symmetries) 


s~i (t j r p\ 

G_{U, 1 ) 


Superconductors with spin-orbital coupling and real pairing 
icjni • <rc;/ + iCjU2 ■ (TCji + ic\nz ■ crc k i + (c^c^,- — cn&fj) (symmetry: time reversal symmetry) 


f ' frp\ 7 

U- [1 ) — Z4 


Superconductors with S z conserving spin-orbital coupling and real pairing 
\c[a z Cj + (c-fiCij — &iiCfj) (symmetry: time reversal and S z spin rotation symmetries) 


G + (U T) — U(l) x Zo 


Superconductors with coplanar spin order and real pairing cjni ■ ac\ + c^n 2 ■ a&j + (c-fiCij — c-iiC^j) 
(symmetry: a combined time reversal and 180° spin rotation symmetry) 


G+(T) = Zl x Zl 


Superconductors with real pairing and collinear spin order c\a z Cj + (c-fiCij — c^Cfj) 
(symmetry: S z spin rotation and a combined time reversal and 180° S y spin rotation symmetry) 


G;(f/,r) = t/(i)x2' 2 T 


Insulators with coplanar spin order c\n\ ■ <TCi + cjn 2 • crcj 
(symmetry: charge conservation and a combined time reversal and 180° spin rotation symmetries) 


G;(!/ 1 r) = [/(i)»z 2 T 


Superconductors with real triplet S z = paring c\i&ij + ciiC-fj 
(symmetry: a combined 180° S 1 ,, spin rotation and time reversal symmetry, a combined 180° S 1 ,, 

\i7C y it j j i y 

spin rotation and charge rotation symmetry, and S z -spin rotation symmetry) 


G77(U,T,C) 

" 4 — j- V 5 ' J 


Superconductors with time reversal, 
180° ^-spin rotation, and SVspin rotation symmetries 


Gtt(U,T, C) 


Superconductors with real singlet pairing c^iCn — c^idfj 
(symmetry: time reversal and SU(2) spin rotation symmetries) 


G[SU(2),T] 


Superconductors with 180° S'y-spin rotation and S^-spin rotation symmetries 


G-(U,C) 


Superconductors with complex singlet pairing e l)>ij (c^c^,- — cj.iCfj) 
(symmetry: 5(7(2)-spin rotation symmetry) 


SU(2) 


Insulators with spin-orbital coupling and inter-sublattice hopping 
ic\ A m ■ crc lB + ic f jA n 2 ■ <rcj B + ic^ns • crc kB 
(symmetry: charge conservation, time reversal and charge conjugation symmetries) 


Gzt(U,T,C) 



TABLE I: (Color online) Electron systems and their full symmetry groups Gf. The groups are defined in Table IV. The 
symmetry group symbols have the following meaning: for example, G„ (U, T, C) is a symmetry group generated by N (the 
U(l) fermion number conservation or spin rotation), T (the time reversal), and C (the charge conjugation or 180° spin rotation). 
The ± subscripts/superscripts describe the relations between the transformations N, T, and/or C (see Table IV). Some times, 
when we describe the symmetry of a fermion system, we do not include the fermion number parity transformation (— ) JV in the 
symmetry group G. Here the full symmetry group G/ does include the fermion number parity transformation (— ) JV . So the full 
symmetry group of a fermion system with no symmetry is G/ = Z^ generated by the fermion number parity transformation. 
Gf is a i?2 extension of G: G = Gf/Z(. Free electron systems with symmetry Gll((7,T,G) actually have a higher symmetry 
G[SU(2), T]. Similarly, free electron systems with symmetry G-(U, C) actually have a higher symmetry 5(7(2). The groups in 
blue are Abelian and electron operators form their ID representations. 



(ie stable even without any symmetry), such as topolog- 
ical superconductors with no symmetry. Just like the 
interacting topological ordered phases, the topological 
phases for free fermions are also characterized by their 
gapless boundary excitations. The boundary excitations 
play a key role in the theory and experiments of free 
fermion topological phases. 

For non-interacting fermion systems with time rever- 
sal (generated by T), charge conjugation (generated by 
C), and/or £7(1) (generated by N) symmetries, the total 
symmetry group may not simply be x x £7(1). 
The group can take different forms, depending on the 
different relations between those symmetry operations, 
such as TNT- 1 = N or TNT- 1 = —N. As a result, 
the gapped phases of those fermion systems with differ- 
ent symmetry groups have different classifications. In 
this paper, we use Kitaev's K-theory approach to clas- 



sify the gapped free fermion phases for those different 
symmetry groups. In Table I, we list some electron sys- 
tems and their full symmetry group Gf. In Table II 
and Table III the 10 classes 25 ' 28 of gapped free fermion 
phases protected by those many-body symmetry groups 
(and many other symmetry groups) are listed. Here we 
have assumed that the fermions form one irreducible rep- 
resentation of the full symmetry group. The result will 
differ if the fermions contain several distinct irreducible 
representations of the full symmetry group (see section 
HIE). In Ref. 25,28, the 10 classes of gapped free fermion 
phases are already associated with many different many- 
body symmetries of electron systems. In this paper, we 
generalize the results in Ref. 25,28 to more symmetry 
groups. 

We note that electron systems, with TNT- 1 = N, 
only realize a subset of the possible symmetry groups. 
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TABLE II: (Color online) Classification of the gapped phases of non-interacting fermions in ci-dimensional space, for some 
symmetries. The space of the gapped states is given by C p +d mod 2, where p depends on the symmetry group. The distinct 
phases are given by 7ro(C p +d mo( j 2). "0" means that only trivial phase exist. Z means that non-trivial phases are labeled by 
non-zero integers and the trivial phase is labeled by 0. The groups in red can be realized by electron systems (see Table I). 
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TABLE III: (Color online) Classification of gapped phases of non-interacting fermions in d spatial dimensions, for some 
symmetries. The space of the gapped states is given by R p ~d mod 8, where p depends on the symmetry group. The phases are 
classified by iro(R p -d mod s)- Here Z2 means that there is one non-trivial phase and one trivial phase labeled by 1 and 0. The 
groups in red can be realized by electron systems (see Table I). 



The emergent fermion (such as the spinon in spin liquid) 
may realize other possible symmetry groups, since their 
symmetries are described by projective symmetry groups 
(PSG) which can be different for different topologically 
ordered states. 41 

The p = line in Table II classifies two types of 
electron systems: (1) Insulators with only fermion num- 
ber conservation (which includes integer quantum Hall 
states). (2) Superconductors with only S z spin rotation 
symmetry which can be realized by superconductors with 
collinear spin order. The p = 1 line in Table II classifies 
superconductors with only time reversal and S z spin ro- 
tation symmetry (full symmetry group G_(U,T)), which 
can be realized by superconductors with real pairing and 
S z conserving spin-orbital coupling. 

In Table III, the p = column classifies electronic in- 



sulators with coplanar spin order (full symmetry group 
G+(U,T) which contains the charge conservation and 
a time reversal symmetry). The p = 1 column classi- 
fies electronic superconductors with coplanar spin order 
and real pairing (full symmetry group G+(T) which con- 
tains a time reversal symmetry). The p = 2 column 
classifies electronic superconductors with non-coplanar 
spin order (full symmetry group "none"). The p = 3 
column classifies electronic superconductors with spin- 
orbital coupling and real pairing (full symmetry group 
G_ (T) which contains the time reversal symmetry) . The 
p = 4 column classifies electronic insulators with spin- 
orbital coupling (full symmetry group GZ (U, T) which 
contains the charge conservation and the time reversal 
symmetry). The p = 5 column classifies electronic in- 
sulators on bipartite lattices with spin-orbital coupling 
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and only inter-sublattice hopping (full symmetry group 
GZ^(U,T,C) which contains the charge conservation, 
the time reversal symmetry, and a charge conjugation 
symmetry). The p — 6 column classifies electronic spin 
singlet superconductors with complex pairing (full sym- 
metry group SU(2)). The p = 7 column classifies elec- 
tronic spin-singlet superconductors with real pairing (full 
symmetry group G[SU(2),T] which contains the SU(2) 
spin rotation and the time reversal symmetry). 

In this paper, we will first discuss a simpler case where 
fermion systems have only the U(l) symmetry. Then we 
will discuss a more complicated case where fermion sys- 
tems can have time reversal, charge conjugation, and/or 
U(l) symmetries. The classification of the gapped phases 
with translation symmetry and the classification of non- 
trivial defects with protected gapless excitations will also 
be studied. 



under the unitary transformation l^ lxl ® | with 

\ W mxm J 

Vlxl <= U(l) and W mxm G U(m). Thus the space C 

of the hermitian matrix satisfying H 2 = 1 is given by 

Um + tn)/U{l) x U(m), which, in n — > oo limit, has 

a form 



C = 



m) 



U(l) x U(m] 



x Z 



(3) 



Clearly 7r (Co) = Z, which recovers the result obtained 
above using a simple argument: the O-dimcnsion gapped 
phases of free conserved fermions are labeled by integers 



B. The properties of classifying spaces 



II. GAPPED FREE FERMION PHASES - THE 
COMPLEX CLASSES 

A. The d = case 

Let us first consider a O-dimension free fermion sys- 
tem with 1 orbital. How many different gapped phases 
do we have for such a system? The answer is 2. The 
2 different gapped phases are labeled by m = 0,1: the 
m = gapped phase correspond to the empty orbital, 
while the m = 1 gapped phase correspond to the occu- 
pied orbital. But 2 is not the complete answer. We can 
alway add occupied and empty orbitals to the system and 
still regard the extended system as in the same gapped 
phase. So we should consider a system with n orbitals 
in n — > oo limit. In this case, the O-dimension gapped 
phases are labeled by an integer m in Z, where m (with 
a possible constant shift) still corresponds to the number 
of occupied orbitals. 

Now let us obtain the above result using a fancier 
mathematical set up. The single-body Hamiltonian of 
the n-orbital system is given bynxn hermitian matrix 
H. If the orbitals below a certain energy are filled, we 
can deform the energies of those orbitals to —1 and de- 
form the energies of other orbitals to +1 without closing 
the energy gap. So, without losing the generality, we can 
assume the H to satisfy 



H 2 



1 



Such a hermitian matrix has a form 



H = U n 



Hxl 




Ul 



(1) 



(2) 



where n = I + m and U nxn inU (n) is an n x n unitary 
matrix. But U nxn is not an one-to-one labeling of the 
hermitian matrix satisfying H 2 = 1. To obtain an onc-to- 



one labeling, we note that 







is invariant 



The space Co is the complex Grassmannian - the space 
formed by the subspaces of (infinity dimensional) com- 
plex vector space. It is also the space of the hermitian 
matrix satisfying H 2 = 1. Actually, Co is a part of a 
sequence. More generally, a space C p can be defined by 
first picking p fixed hermitian matrices 7j, i — 1, 2, ...,p, 
satisfying 

HI j = 2<%- (4) 

Then C p is the space of the hermitian matrix satisfying 

H 2 = l, li H = -H li , i = l,..., p. (5) 

To see what is the G\ space, let us choose 71 — Inxn 
that satisfy 7^ = 1. But for such a choice, we cannot 
find any H that satisfy 71 H = —Hji. Actually, we have 
a condition on the choice of 71. We must choose a 71 
such that 71 if = —Hji and H 2 = 1 has a solution. 
So we should choose 71 = <r x <£> I nxn - We note that 71 
is invariant under the following unitary transformations 
e ia*®A nx „ e i<7°®B nxn g jj^ x u( y ) ( whcrc Anxn and 

B nxn are hermitian matrices). Then H satisfying H 2 = 1 
and 71 if + f?7i = has a form 

H= (6) 

e icr x (g) J 4„ >< „ e icr (g)B„ x „ / nxn )e -i<T0 ® B " x " e -io ' X ®" 4 " x " 

whose positive and negative eigenvalues are paired. We 
see that the space C\ is U(n) x U(n)/U(n) = U(n). 
To construct the C 2 space, we can choose 71 = a x ® 



and 72 = a v 



Then H satisfying H 2 = 1 



and 7,ii + H 7$ = 0, i = 1, 2, has a form 



H = a° ® U nxn 




(7) 

where n = I + m and C/„ x „ G C/(n). We see that the 
space C2 = Co- 
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To construct the C3 space, we can choose 71 = a x 
Inxn, 72 = a v ® 4xni an d 73 = c z ® J BXn . But for such 
a choice, the equations 73 i? + if 7, = 0, i = 1,2,3, and 
H 2 = 1 has no solution for H . So we need to impose the 
following condition on Tj's: 

The equations 7,-!? + i?7i = 0, H 2 = 1 

have a solution for iJ. (8) 

(Later, we will see that such an condition has an amazing 
geometric origin.) Let us choose 71 = a x <g> a x <g> I nxn , 
^2 = <J V ® & x ® I n xn: an d 73 = <t z ® cr 21 ® I nX n instead. 
Then 77 satisfying # 2 = 1 and j.H + H^ = 0, i = 1, 2, 3, 
has a form 



matrix _ff , whose continuous limit has a form: 



i? = e 1 



i(T u <g>CT u <g>.B,, x „ / o _z 



(cr U (8)(7 Z ®/„xn)X 



(9) 



We find that C 3 = d. 

Now, it is not hard to sec that C p = C p +2 which leads 
to ir d (C p ) = TTd(C p+2 ). Thus 



Z, p = mod 2, 
{0}, p = 1 mod 2. 



(10) 



C. The d / cases 

Next we consider a d-dimension free conserved fermion 
systems and their gapped ground states. Note that the 
only symmetry that we have is the U(l) symmetry asso- 
ciated with the fermion number conservation. We do not 
have translation symmetry and other symmetries. 

To be more precise, our e?-dimension space is a ball 
with no non-trivial topology. Since the systems have a 
boundary, here we can only require that the "bulk" gap 
of the fermion systems are non-zero. The free fermion 
system may have protected gapless excited at the bound- 
ary. (Requiring the fermion systems to be even gapped 
at boundary will only give us trivial gapped phases.) We 
will call the free fermion systems that are gapped only 
inside of the d dimensional ball as "bulk" gapped fermion 
systems. A "bulk" gapped fermion system may or may 
not be gapped at the boundary. 

Kitacv has shown that the space of such "bulk" 
gapped free fermion systems is homotopically equivalent 
of the space of mass matrices of a d-dimensional 
Dirac equation: n n (C^) = TTn(C^). 25 In the following, 
we will give a hand- waving explanation of the result. 

To start, let us first assume that the fermion system 
has the translation symmetry and the charge conjuga- 
tion symmetry. We also assume that its energy bands 
have some Dirac points at zero energy and there are no 
other zero energy states in the Brillouin zone. So if we 
fill the negative energy bands, the single-body gapless 
excitations in the system are described by the hermitian 



H = ^7ti<9i 



(11) 



where we have folded all the Dirac points to the k = 
point. Without losing generality, we have also assumed 
that all the Dirac points have the same velocity. Since 
idi is hermitian, thus ji, i = l,...,d are the hermitian 
7-matrices (of infinity dimension) that satisfy eqn. (4). 

When d = 1, do we have a system that has 71 = Inxn? 
The answer is no. Such a system will have n right- 
moving chiral modes that cannot be realized by any pure 
1-dimcnsional systems with short-ranged hopping. In 

'hxl 
l Imx 

I = m (the same number of right-moving and left-moving 
modes). So the allowed 7! always satisfy the condition 
that H 2 = 1 and 71 H + H 71 = has a solution for H . 
We see that the extra condition eqn. (8) on 7^ has a very 
physical meaning. 

Now we add perturbations that may break the transla- 
tion symmetry. We like to know how many different ways 
are there to gap the Dirac point. The Dirac points can 
be fully gapped by the "mass" matrix M that satisfies 



fact 71 must have a form ji — 



with 



TiM + Mji = 0, M f = M. 



(12) 



To fully gap the Dirac points, M must have no zero eigen- 
values. Without losing generality, we may also assume 
that M 2 — 1. (Since the Dirac points may have dif- 
ferent crystal momenta before the folding to k = 0, we 
need perturbations that break the translation symmetry 
to generate a generic mass matrix that may mix Dirac 
points at different crystal momenta.) The space of 
such mass matrices is nothing but Ca introduced before: 
Cf = C d .^ 

So the different ways to gap the Dirac point form a 
space Cd- The different disconnected components of C d 
represent different gapped phases of the free fermions. 
Thus, the gapped phases of the conserved free fermions 
in ^-dimensions are classified by TTo(Cd), which is Z for 
even d and for odd d. The non-trivial phases at d = 2 
are labeled by Z, which are the integer quantum Hall 
states. The results is summarized in Table II. 



III. GAPPED FREE FERMION PHASES 
REAL CLASSES 



THE 



When fermion number is not conserved and/or when 
there is a time reversal symmetry, the gapped phases of 
non-interacting fermions are classified differently. How- 
ever, using the idea and the approaches similar to the 
above discussion, we can also obtain a classification. In- 
stead of considering hermitian matrices that satisfy cer- 
tain conditions, we just need to consider real antisym- 
metric matrices that satisfy certain conditions. 
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A. The d — case — the symmetry groups 

Again, we will start with d — dimension. In this 
case, a free fermion system with n orbitals is described 
by the following quadratic Hamiltonian 

H = ^H lj c\c J + ^2 l [G lJ c l c j + h.c], i,j = 1, ...,n 

ij ij 

(13) 

Introducing Majorana fermion operator fji, I = 1, 2n, 
{fji,fjj} = 25u, 4 = Vi, (14) 
to express the complex fermion operator df. 
1 



Ci = -{fhi + 



we can rewrite H as 



(15) 



(16) 



ij 



where A is a real antisymmetric matrix. For example, 
for a 1-orbital Hamiltonian H = e{&c — |), we get A = 

C -.1 

If the fermion number is conserved, H commutes with 
the fermion number operator 



IJ 



where 



Q = e®J, Q 2 = -l, s = -ia y . (18) 

[H, N] =0 requires that 

[A,Q]=0. (19) 

Such matrix A has a form A = a ® i? a + e <£> H s , where 
i7 s is symmetric and if a antisymmetric. We can convert 
such an antisymmetric matrix A into a hermitian ma- 
trix H = H s + iH a , and reduce the problem to the one 
discussed before. 

The time reversal transformation T is antiunitary: 
TiT^ 1 = — i. Since T does not change the fermion num- 
bers, thus TdiT^ 1 — Uijdj, where U is an unitary matrix. 
In terms of the Majorana fermions, we have 

ffoiT^ 1 = RcUij^j - ImUijfj 2 j+i, 
f^+if- 1 = -ReU ijV2j +i - hnUijihj. (20) 

Therefore, we have 

ffiif- 1 =T ij fij, T = a 3 (giReU - a 1 <g> ImZ7. (21) 



We see that in the Majorana fermion basis, U — > a 3 ® 
ReU — a 1 <g> ImU = T and i -> e ® J. We indeed have 
T(e® J) = -T(e® J). 

For fermion systems, we may have T 2 = St = ±- 
In fact s T = — for electron systems. This implies that 
T 2 CiT~ 2 = STCi and T 2 = sy. The time reversal invari- 
ance fkf- 1 = H implies that T T AT = -A, where T T 
is the transpose of T. We can show that 

T T T = (a 3 <g> Re[/ f - ct 1 <g> ImC/ t )(c7 3 <g> Ret/ + cr 1 <g> Imf/) 
= cr° <g) (ReL^Ret/ - ImZ/tlmt/) 
- e <g> (Ref/^mC/ + Im^Ref/) 
= o- ® 7, (22) 

where we have used ReWReU — ImWlmU = I and 
Re?7 t Im[/ + Im£/ t ReJ7 = for unitary matrix U. There- 
fore T T = T _1 and 



AT = -TA, T 2 = st- 



(23) 



Also, for fermion systems, the time reversal transforma- 
tion T and the U(l) transformation TV may have a non- 
trivial relation: TV^T" 1 = e s ^ iej \ 



TNT- 1 = -s UT N. This gives us 
TQ = s UT QT. 



SUT 



±, or 



(24) 



The charge conjugation transformation C is unitary. 
Since C changes Ci to c\, thus CciC^ 1 = UijCj, where U 
is an unitary matrix. In terms of the Majorana fermions, 
we have 

Cr)nC~ x = ReUijf)2j + ImUijfcj+i, 
Cfj 2 i+iC~ l = -ReUijfi2j+i + ImUijfcj. (25) 

Therefore, we have 

CfitC- 1 = djfjj , C = a 3 ® ReU + a 1 <g> ImU. (26) 

Again, we can show that C T = C _1 . 

For fermion systems, we may have C 2 — , sc = ±, 
which implies that C 2 CiC~ 2 = sgCi and C 2 = so The 
charge conjugation invariance CHC^ 1 = H implies that 
A satisfies 



CA^CA, C 2 = s c . 
Since CiVC" 1 = - A, we have 

CQ = -QC. 



(27) 



(28) 



However, the commutation relation between f and C has 
two choices: TC = s^ c CT 7 stc — ±, we have 



(29) 



We see that when we say a system has C/(l), time re- 
versal, and/or charge conjugation symmetries, we still 
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Symmetry groups 


Relations 


U(l), SU(2) 




Gs C (C) 


C = s c , s c =±. 


Gs C (U,C) 


C — s c , C e C = e , sc = ±- 


Gs T (T) 


f 2 = 4, «T = ±. 


G^ T (C/,T) 


f e ieNf-l = e s UT iBN^ f 2 = S^, S UT ,S T =±. 




f 2 =s§, C 2 = sg, Cf=(s T r c )fC, s TC ,s T ,s c = ±. 


GZTc TC (U,T, C) 


T 2 = 4, ^=4, 6f={4c)f&, ST,SC,SUT,STC=±- 



TABLE IV: Different relations between symmetry transformations gives rise to 36 different groups that contain U(l) (repre- 
sented by 17), time reversal T, and/or charge conjugation C symmetries. 



do not know what is the actual symmetry group of the 
system, since those symmetry operations may have differ- 
ent relations as described by the signs St, sc: s ut, s tco 
which lead to different full symmetry groups. Because 
symmetry plays a key role in our classification, we cannot 
obtain a classification without specifying the symmetry 
groups. We have discussed the possible relations among 
various symmetry operations. In the Table IV, we list 
the corresponding symmetry groups. 

We like to point that that some times, when we de- 
scribe the symmetry of a fermion system, we do not in- 
clude the fermion number parity transformation (— ) N in 
the symmetry group G. However, in this paper, we will 
use the full symmetry group Gf to describe the symme- 
try of a fermion system. The full symmetry group Gf 
does include the fermion number parity transformation 
(— ) N . So the full symmetry group of a fermion system 
with no symmetry is Gf — z£ generated by the fermion 
number parity transformation. Gf is actually a Z 2 ex- 
tension of G: G = Gfjz{. It is a projective symmetry 
group discussed in Ref. 41. 

In the following, we will study the symmetries of var- 
ious electron systems, to see which symmetry groups 
listed in Table IV can be realized by electron systems. 

For insulators with non-coplanar spin order SH = 
c\ni ■ crci + cjri2 • crcj + cj.n.3 • crdk, the full symmetry 
group is Gf — U{\) generated by the total charge Nc- 

For superconductors with non-coplanar spin order 
SH = clni-aci+cjn2-(TCj+cln3-ack + (c^iCij-ciiCfj), 
the full symmetry group is reduced to Gf = Z% generated 
by the fermion number parity operator Pf = (—) Nc . We 
note that the full symmetry group of any fermion system 
contain Z^ as a subgroup. So we usually, use the group 
Gf/Z^ to describe the symmetry of fermion system, and 
we say there is no symmetry for superconductors with 
non-coplanar spin order. But in this paper, we will use 
the full symmetry group G f to describe the symmetry of 
fermion systems. 

For insulators with spin-orbital coupling SH = icjni • 
ercV + icjn 2 • ercy + ic^r^ • crdk', they have the charge 



conservation Nc and the time reversal Tph y symmetries. 
The time reversal symmetry is defined by 

TphyC Q ,iT ph y = e Q( 3C/3 : i, TphycJ^Tpjjy = e a pc.p^. (30) 

We can show that 

TphycjerCjTphy = -cJctCj, 

T phy N c f-^ = N c , T p 2 hy = R^ (31) 

Thus SH is invariant under Tph y and e l9Nc . Let T — 
Tp hy and N — N c , we find that 

= e -itW f2 = ( _)JV (32) 

which define the full symmetry group GZ(U,T) of an 
electron insulator with spin-orbital coupling. 

For superconductors with spin-orbital coupling and 
real pairing SH = icfni ■ ercv + ic]«2 ■ trcji + icj.n 3 • 
ercv + (c-fiCij — ciidfj), they have the time reversal sym- 
metry Tphy. Setting T = T p h y and N = Nc, we find 

f 2 = (33) 

which defines the full symmetry group G_(U,T) = Z4 
of superconductors with spin-orbital coupling and real 
pairing. 

For superconductors with S z conserving spin-orbital 
coupling and real pairing SH = ic\a z Cj + (&|-jCj.j — c^dfj) , 
they have the time reversal T p h y and S z spin rotation 
symmetries. Setting T = T ph y and N = 2S Z , we find 

f e i6Nf-i = e ifw f2 = ( _yv ; (34) 

which defines the full symmetry group G~Z(U,T) = 
U(l) x Z2 of superconductors with S z conserving spin- 
orbital coupling and real pairing. 

For superconductors with real pairing and coplanar 
spin order SH — cjni ■ ercj + cjn 2 ' + (cfjCjj — cjjCfj), 
they have a combined time reversal and 180° spin ro- 
tation symmetry. The spin rotation is generated by 

<Sa = J2i ^ c \ a ' ac i' a = x iV, z - We have 

Tphy S a T h = —S a . (35) 
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The Hamiltonian SH is invariant under T — e 17rS f p h y . 
Since T 2 = 1, the full symmetry group of superconduc- 
tors real pairing and coplanar spin order is G+(T) — 
Z 2 x Z{. 

For superconductors with real pairing and collinear 
spin order SH — c\a z Cj + (cfiCjj ~ c^iCfj), they have 
the S z spin rotation and a combined time reversal and 
180° S y spin rotation symmetry. The Hamiltonian SH 

is invariant under T = e 17rSV T p hy and S z spin rotation 
N = 2S Z . We find that 



Te ieN T 



-i8N 



(-) 



A'* 



(36) 



which define the full symmetry group G + (U, T) of super- 
conductors with real pairing and collinear spin order. 
For insulators with coplanar spin order cjni • erci + 

CjTi2 ■ crcj, they have the charge conservation and a com- 
bined time reversal and 180° spin rotation symmetries. 
The Hamiltonian SH is invariant under T = e'^Jphy 
and the charge rotation N = Nc- We find that 



fe ieN f 



-i8N 



rp'l 



(") 



N 



(37) 



which define the full symmetry group G7(U,T) of insu- 
lators with coplanar spin order. 

For superconductors with real triplet S z = paring 
SH = c^i&ij + c^ic^j, they have a combined time reversal 
and charge rotation symmetry a combined 180° S y spin 
rotation and charge rotation symmetry, and the 5 z -spin 
rotation symmetry. The Hamiltonian SH is invariant un- 
der f = e'f^fphy, C = e 1 f^cr e i^ J)j and the Sz spin 

rotation N = 2S Z . We find that 



the Hamiltonian is invariant under T = T p hy, C = c 
and N = 2S Z . We find that 



Ce^C' 1 = e- lHiV , Te T 

f 2 = C 2 = cf = fc. 



8Nrfi-l _ i8N 



(40) 



which define the full symmetry group G^(U,T,C) of 
superconductors with the time reversal, the 180° S^-spin 
rotation, and the S^-spin rotation symmetries. For free 
electrons with the 180° S'y-spin rotation, and the S^-spin 
rotation symmetries, they actually have the full SU(2) 
spin rotation symmetry. So the above systems are also 
superconductors with real pairing and the SU(2) spin 
rotation symmetry. Similarly, for superconductors with 
complex pairing and the SU(2) spin rotation symmetry, 
the symmetry group is SU(2), or G-(U,C). 

For insulators with spin-orbital coupling and only 
inter-sublattice hopping H = icj^ni • crc iB + ic^n 2 • 

crcj B + icj.^n3 • <TCfe B) they have the charge conservation, 
the time reversal and a deformed charge conjugation sym- 
metries. The charge conjugation transformation C* p hy is 
defined as 

C'phyCa.iC'phy — £ Q / 3c] 3jl , A>hyCa,i^phy = e apCp,i- (41) 

We find that 



A>hyCj0'CjC' phy — C^CTCj, CphyTphy — TphyC'phy, 

C phy N c C-' = -N c , C 2 phy = (-fo 



(42) 



The above Hamiltonian is invariant under T = Tphy, N — 

Nc, and C — (— ) NB C p \ ly , where Nb is the number of 
electrons on the B-sublatticc. We find 



i8N 



fjOfif-1 = e i8N^ 

f 2 = 1, C 2 = 1, cf = (-ffc 



(38) 



Ce^C- 1 = e~ ieN , 

f e i8Nf-l = c -i8N^ 

f 2 = (-)*, C 2 = (-)*, cf = fc 



(43) 



which define the full symmetry group G++(U,T,C) of 
superconductors with real triplet S z — paring. 

For superconductors with real triplet S z = paring 
and collinear spin order SH — CiO z Ci + (&|-jCy + c^c^), 
they have a combined time reversal and 180° S v spin 
rotation symmetry, and the 5 z -spin rotation symmetry. 
The Hamiltonian SH is invariant under T = e 1 ^ Sy T p i iy , 
and the S z spin rotation N — 2S Z 



We find that 



-iON 



(39) 



which define the full symmetry group G^_(U, T) of super- 
conductors with real triplet S z — paring and collinear 
spin order. 

For superconductors with the time reversal, the 180° 
Sy-spin rotation, and the S^-spin rotation symmetries, 



which define the full symmetry group G~_(U,T,C) of 
insulators with spin-orbital coupling and only inter- 
sublattice hopping. The above results for electron sys- 
tems and their full symmetry groups Gf are summarized 
in Table I. 



B. The d = case — the classifying spaces 

The hermitian matrix iA describes single-body excita- 
tions above the free fermion ground state. We note that 
the eigenvalues of \A are ±e$. The positive eigenval- 
ues \ei\ correspond to the single-body excitation energies 
above the many-body ground state. The minimal |e$| rep- 
resents the excitation energy gap of H above the ground 
state (the ground state is the lowest energy state of H) . 
So, if we are considering gapped systems, |e,| is alway 
non-zero. We can shift all to ±1 without closing the 
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gap and change the (matter) phase of the state. Thus, 
we can set A 2 = —1. 

In presence of symmetry, A should also satisfy some 
other conditions. The space formed by all those A's is 
called the classifying space. Clearly, the classifying space 
is determined by the full symmetry group G/. In this 
section, we will calculate the classifying spaces for some 
simple groups. 

If there is no symmetry, then the real antisymmetric 
matrix A satisfies 

A 2 = -1. (44) 

The space of those matrices is denoted as Rq, which is 
the classifying space for trivial symmetry group. 

If there is only the charge conjugation symmetry (full 
symmetry group = G 8C (C)), then the real antisymmetric 
matrix A satisfies 

A 2 = -l, AC = CA, C 2 = s C - (45) 

For sq = +j since C commute with A and C is symmet- 
ric, we can always restrict ourselves in an eigenspace of 
C and C can be dropped. Thus the space of the matrices 
is Rq, the same as before. For sc = — , we can assume 
C = e<E)I- In this case, A has a form A = a°®H a +e<E)H s , 
where H s = Hj and H a = — H£ . Thus, we can convert 
A into a hermitian matrix H = H s + iH a , and the space 
of the matrices is Co- 

If there are ?7(1) and charge conjugation symmetries 
(full symmetry group = G SC (U,C)), then the real anti- 
symmetric matrix A satisfies 

A 2 = -1, AQ = QA, AC = CA, QC = -CQ, 
Q 2 = -l, C 2 = s c - (46) 

For sq = +, we can assume C — g z ® I and Q = e ® I. 
Since Q and C commute with A, we find that A must 
have a form A = er° <E) A, with A 2 = — 1. Thus the space 
of the matrices A, and hence A, is Rq. 

For sc = — , we can assume C — e <£> a x ® I and 
Q = e <£> a z <£> I. We find that A must have a form 
A = a ®<T ®Ho + e®cr () ®H 1 +<j z ®e®H 2 +<j x ®e®H 3 , 
where Hq = —Hq and Hi = Hf, i = 1,2,3. Now, we 
can view a ® cr° as 1, e <g> a as i, a z <£> e as j, and 
a x (g> e as k. We find that i,j,k satisfy the quaternion 
algebra. Thus A can be mapped into a quaternion matrix 
H = H + \Hi + }H 2 + kiJ 3 satisfying W = -H and 
H 2 = —1. The quaternion matrices that satisfy the above 
two conditions has a form 

H=e x —\I nxn e- x — (47) 

where Xj ixn = —X nxn is a quaternion matrix. e x " x " 
form the group Sp(n). However, the transformations 
c -4„ x „ + iB„ x „ k CC p S H nxn unchanged, where A nxn is a 
real antisymmetric matrix and B nxn is a real symmetric 
matrix. e" 4 " x " +lS " x ™ form the group U(n). Thus, the 
space of the quaternion matrices that satisfy the above 



two conditions is given by Sp(n)/U(n). Such a space is 
the space Rq which will be introduced later. 

When sc = — 1, we can view Q as the generator of S z 
spin rotation, and C as the generator of S x spin rotation 
acting on spin- 1/2 fermions. In fact e *® and e 9xC in this 
case generate the full SU(2) group. So when = the 
free spin- 1/2 fermions with U(l) x Z^ 7 symmetry actually 
have the full SU(2) spin rotation symmetry. Therefore 
G-(U,C)~SU(2). 

If there is only the time reversal symmetry (full sym- 
metry group = G ST (T)), then A satisfies 

A 2 = -1, A Pl + Pl A = 0, pj = s T , Pi = T. 

(48) 

The space of those matrices is denoted as R^ for st = — 1 
and i?J for St = 1- 

If there are time reversal and U(l) symmetries (full 
symmetry group = G S S ^, T (U, T)), then for sjjt = — , A 
satisfies 

A 2 = -1, Api + piA = 0, p\ = s T , p\ = s T , 

Pl = T, P2 =TQ. (49) 

The space of those matrices is denoted as i?° for st = + 
and Rq for st = —■ 

For sjjt — +, Q commute with both A and T. Since 
Q 2 = — 1, we can treat Q as the imaginary number i and 
convert both A and T to complex matrices. To see this, 
let us choose a basis in which Q has a form Q = e ® /. 
In this basis A and T become A = <r° ® A 2 + e <£> A\ 
and T = <r ® Ti + e (£> T2 , where is symmetric and 
A 2 is antisymmetric. Let us introduce complex matrices 
H = -Ai + \A 2 and f = T x + \T 2 for s T = + or f = 
-T 2 + ill for s T = -. From A 2 = -1, T 2 = s T , and 
AT = -TA, we find 

ff 2 = l, HT + TH = 0, f 2 = l. (50) 

Also A T = — A allows us to show W = H. For a fixed T, 
the space formed by -ff's that satisfy the above conditions 
is G\ introduced before. This allows us to show that the 
space of the corresponding matrices A is G\ for s T = ±, 

SUT = +■ 

If there are time reversal and charge conjugation sym- 
metries (full symmetry group = G S S ^ C (T,C)), then for 
stc = —■, A satisfies 

A 2 = -1, Api + piA = 0, p\ = s T , p\ = -s T sc, 
Pl =T, P2 =TC. (51) 

The space of the matrices A is R\ for st — +, sc — +; 
i?2 for s t = +, sc = — ! #1 f° r s t = —, sc = +; and Rq 
for s T — — , s c = — . For s TC = +, C will commute with 
both ^4 and T. We find space of the matrices A to be i?5 
for = +, sc = +; i?o for st = — , sc = +; and Ci for 

s T = ±, sc = -■ 

If there are U(l), time reversal, and charge conjugation 
symmetries (full symmetry group = G S S ^^, TC (U, T, C)), 
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then for stc — s ut = A satisfies 

A 2 = -1, Ap, + p,A = 0, p\= p\ = s t , pi = -s T s c , 
Pi=T, p 2 =TQ 1 P3 = TC. (52) 

The space of the matrices A is R\ for st — +, sc = +; 
i?3 for st = +, s c = — ; i?i for st — — , s c = +; and i?g 
for s T = sc = -■ 

For sjtt = — , stc = +, A satisfies 

A 2 = — 1, Ap 4 + p 4 v4 = 0, p\= p\ = st, p\ = stSc, 
Pi=T, p 2 =TQ, p 3 =TQC. (53) 

The space of the matrices A is R\ for st = +, sc = +; 
i?3 for st = +, sc = ~ ; R\ for st = — , sc — +; and i?g 
for st = — j sc = — • 

For syT = +j stc = ~, ^4 satisfies 

A 2 = -1, Api + p t A = 0, p\ = s T , p\ = p\ = -s T s c , 
p x =T, P2 = TC, p 3 = TCQ. (54) 

The space of the matrices A is R 2 for st = +, sc = +\ 
i?3 for st = +, sc = — ; R\ for st = — , sc = +; and Rq 
for st = -, sc = -■ 

For S[/t = +, stc = +, we nncl that A satisfy 

A 2 = -1, Ap 4 + p 4 v4 = 0, p 2 = -st, p\ = pi = s T sc, 
Pi = TQ, P2 = TC, P3 = TCQ. (55) 

We see that the matrices A form a space R\ for st = 
+, sc = +; Rq for st = +, sc = -; -R 2 for st = 
-, s c = +; and R° for s T — — , s c = — ; 



We can check that such pi satisfy the following Clifford 
algebra Cl(p+l,q + 2) 



PjPi + PiPj 



Pi 



1. 



0, 



,p+l 



p\ = 



- 1. 



i=p+2,...,p+q+3 



(60) 



If we fix i 7^ p + q + 2, then the space formed by 
v4 = p p+q+ 2 satisfying the above condition is given by 
RpXi- The above construction gives rise to a map from 
Rp" —} Rp%x- Since A — p p + q +2 satisfying eqn. (60) must 
has a form A®a z , with A satisfying eqn. (58). This gives 
us a map R q p %\ -> R*. Thus R q p %\ = R q . 

We can also consider real symmetric matrix A that 
satisfies (for fixed real matrices pi, i = 1, ...,p) 



A = p p +i, pjPi + piPj 



o, Pi 



i=l,...,p+l 



The space of those matrices is denoted as R p . 
In the following, we are going to show that 



i?n — R, 



q+2- 



(61) 



(62) 



From A £ Rq that satisfies the Clifford algebra CZ(0, q - 
1) 



A = 



Pq+1> 



Pi Pi + PiPj = 



0. 



Pi 



,9+1 



1. 



(63) 



we can define 



C. The properties of classifying spaces 

In general, we can consider a real antisymmetric matrix 
A that satisfies (for fixed real matrices pi, i = 1, ...,p + q) 



A = p P+q +i, PjPi + piPj 



0. 



i— 1,... ,p 



i, p} = 



- 1. 



i=p+l,...,p+q+l 



The space of those A matrices is denoted as R q . 
Let us show that 

i\p — n p+ i- 



(56) 



(57) 



From A e R q that satisfies the following Clifford algebra 
Cl(p,q + 1) 



A = Pp+q+l, ft/ft ■ P,P, 



0, 



ft- ■ 

2—1, 

we can define 



i, # 



i=p+l,...,p+q+l 



1. 



(58) 



ft 
Pi 



i=l,...,p 



ft ®CT Z , pp+i = I (g)<7 x , 



i=p+2,...,p+q+2 



Pi^®a z , p p+q+ 3 = I <8>e. (59) 



Pi= Pi<Z>e, p q+2 = /® er z , p q+3 = I <E) a x . 

i=l,...,q+l 

(64) 

we can check that p^ form the Clifford algebra C7(g+3, 0) 



PjPi + PiPj 



0, Pi 



1. 



i=l,...,g+3 



(65) 



If we fix i 5 + 1, then the space formed by A = p 9 +i 
satisfying the above condition is given by R q +2- The 
above construction gives rise to a map from R^ — > Rq+2 ■ 
Since A — p q +\ satisfying eqn. (65) must has a form 
A ® £, with A satisfying eqn. (63). This gives us a map 
R q+2 -> i?g. Thus i2g = i? 9+2 . 

In addition we also have the following periodic rela- 
tions 



D9 _ D9+8 _ r>9 



i?p — Rp+8- 



(66) 



This can be shown by noticing the following 16 dimen- 
sional real symmetric representation of Clifford algebra 
C7(0,8): 



9 3 = e <g) cr z (8 

9 7 = £ (g) £ ® (7° ® (7°, 



°(g)£, 6*2 = £ ® cr Z ® £ ® CT 31 , 
£(g(T z , 6*4 = £(g(7 ;E (g)£(g)Cr , 

cr x (8 e, 6» 6 = £ (g cr 31 (g) cr z (g £, 

^ct ®^ , (67) 



(7 (g IT 
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p mod 8 





1 


2 


3 


4 


5 


6 


7 


R p 


C(i+m) v 1 


O(n) 


0( 2 n) 


U(2n) 


Sp(i + m) v ? 


Sp(n) 


Spin) 


t7(n) 


0(l)xO(m) A ^ 


U(n) 


Spin) 


Sp(0xSp(m) A L 


U(n) 


O(n) 


TTo(-Rp) 


z 


z 2 


z 2 





Z 











Tl(-Rp) 


z 2 


z 2 





z 











z 


vr 2 (i?p) 


z 2 





z 











z 


z 2 


7T 3 (iip) 





z 











z 


z 2 


z 2 


7r 4 (-Rp) 


z 











z 


z 2 


z 2 





^(Rp) 











z 


z 2 


z 2 





z 


ire(Rp) 








z 


z 2 


z 2 





z 










z 


z 2 


z 2 





z 









TABLE V: The spaces R p and their homotopy groups nd(Rp)- 



which satisfy 



0, 



We find that 6 
a anticommute with 6, 
eqn. (58), we can define 



¥i 



i=0,...,S 



From A £ R q that satisfies 



Pi 



Pi = 



i— 1,... ,p 



Pp4 



J® 0i 



i=l,...,8 



i=p+9,...,p+g+9 

We can check that such pi satisfy 



(69) 



Pjpi + piPj 



¥3 



0, Pt 



i=l,...,p+8 



- 1. 



(70) 

=p+9,...,p+g+9 

If we fix i 7^ p + q + 9, then the space formed 
by A = p p+q+ g satisfying the above condition is given 
by Rp +S - The above construction gives rise to a map 
from R q p — > Rp +S - On the other hand, the matrix that 
anticommute with all ft's must be proportional to 0. 
Thus A — Pp+ q +9 satisfying eqn. (70) must has a form 
A Cg) 8, with A satisfying eqn. (58). This gives us a map 
R p+S — > Rp- Thus R p+& = R p - Using a similar ap- 
proach, we can show R p = R p +&. The relation eqn. (57), 
eqn. (62), and eqn. (66) allow us show 



Ftp — Rq-p+2 mod 8- 



So we can study the space R q via the space R q - p +2 mod s- 
Let us construct some of the R p spaces. Rq is formed 
by real symmetric matrices A that satisfy A 2 = 1. Thus 



A has a form O 







0-\ O e 0[l + m). We 



Q(l+m) 



where e a ® L ^^ G 0(n) and e CT *® M ™ x " G 0(n) are the 
transformations that leave p\ unchanged. We see that 
(68) R x = 0(n). The other spaces R p and ir Q (R p ) are listed in 
Table V. Note that for space S'xZ.we have ttq(S x Z) = 
cr° ® 7r (5) x Z. Also O(n) in dividend usually leads to the 
Z 2 in 7r . Otherwise 7r = {0}. 0(n) in dividend can 
give rise to Z 2 because O G O(n) with det(O) = 1 and 
det(O) = —1 cannot be smoothly connected. 0(l + m) in 

o(7)xO(to) ^ oes n0 * ^ ea< ^ to ^ 2 because for O G 0(i + m), 
we can change the sign of dct(O) by multiplying O with 
an element in 0(1) [or 0(mj\. 

For free fermion systems in 0-dimension with no sym- 
metry and no fermion number conservation, the classify- 
ing space is Since 7r (i?Q) = 7ro(-R2) = Z 2 , such free 
fermion systems has two possible gapped phases. One 
phase has even numbers of fermions in the ground state 
and the other phase has odd numbers of fermions in the 
ground state. (Note that the fermion number mod 2 is 
still conserved even without any symmetry.) 

For free electron systems in 0-dimension with time re- 
versal symmetry and electron number conservation (the 
symmetry group GZ(U,T)), the classifying space is Rq. 
Since 7To(i?o) = 7r o(^4) = Z, the possible gapped phases 
are labeled by an integer n. The ground state has 2n 
fermions. The electron number in the ground state is 
always even due to the Kramer degeneracy. 

If we drop the electron number conservation (the sym- 
metry group becomes G_(T)), then the ground state 
will have uncertain but even numbers of electrons. The 
ground state cannot have odd numbers of electrons. This 
implies that free electron systems with only time reversal 
symmetry in 0-dimension has only one possible gapped 
phase. This agrees with 7r (i?o) = 7r (i?3) = {0}, where 
Rq is the classifying space for symmetry group G_(T). 



(71) 



see that R = \J m O(l+m)/0(l)xO(m) - 0(l)x0{m 

Ri is formed by real symmetric matrices A that satisfy 
A 2 — 1 and Api = —p\A with py = a 
has a form 



D. The d ^ cases 



I nxn . Thus A 



A = e°*® M - 



Now let us consider the d ^ cases. Again, let 
us first assume that the fermion system described by 
-<r°®L nXn e -<T*®M nxn H = \ AuViVJ has tlic translation symmetry, as 
well as the time reversal symmetry and fermion number 
, (72) conservation. We also assume that the single-body en- 
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ergy bands of antisymmetric hermitian matrix iA have 
some Dirac points at zero energy and there are no other 
zero energy states in the Brillouin zone. The gapless 
single-body excitations in the system are described by 
the continuum limit of iA: 

d 

iA=iJ2lidi, (73) 

i=l 

where we have folded all the Dirac points to the k = 
point. Without losing generality, we have also assumed 
that all the Dirac points have the same velocity. Since 
di is real and antisymmetric, thus i = l,...,d are real 
symmetric 7-matrices (of infinity dimension) that satisfy: 

7i7j + Ijli = 2<%, 1* = 7». (74) 

Again, the allowed 7$ always satisfy the condition that 
M 2 = -1 and jiM + Mji = has a solution for M. 
Since the time reversal and the U(l) transformations do 
not affect di, therefore the symmetry conditions on A, 
AT + TA = and AQ - QA = 0, become the symmetry 
conditions on the 7-matrices: 

TiT + Tji = 0, 7l Q - Q 7i = 0. (75) 

Now we add perturbations that may break the trans- 
lation, time reversal and the U(l) symmetries, and ask: 
how many different ways are there to gap the Dirac 
points. The Dirac points can be fully gapped by real 
antisymmetric mass matrices M that satisfy 

7i M + Af 7l = 0, (76) 

The resulting single-body Hamiltonian becomes iA — 
iEtib^ + M]. 

If there is no symmetry, we only require the real anti- 
symmetric mass matrix M to be invertible (in addition 
to eqn. (76)). Without losing generality, we can choose 
the mass matrix to also satisfy 

M 2 = -1. (77) 

The space of those mass matrices is given by R° d . 

If there are some symmetries, the real antisymmetric 
mass matrix M also satisfy some additional condition, 
as discussed before: M anticommutes with a set of p + q 
matrices pi that anticommute among themselves with p 
of them square to 1 and q of them square to -1. The 
number of p, q depend on full symmetry group G t . Since 
ji do not break the symmetry, so, just like M, 7 j also 
anticommute pi. So, in total, M anticommutes with a 
set of p + q + d matrices pi and 7^ that anticommute 
among themselves with p + d of them square to 1 and q 
of them square to -1. Those mass matrices from a space 
R q 

n p+d- 

The different disconnected components of Rp +d rep- 
resent different "bulk" gapped phases of the free 
fermions. Thus, the "bulk" gapped phases of the free 
fermions in (^-dimensions are classified by ttq = 
-K {R q - p - d+2 mod a), with (p,q) depending on the sym- 
metry. The results are summarized in Table III. 



E. A general discussion 

Now, let us give a gneral discussion of the classifying 
problem of free fermion systems. To classify the gapped 
phases of free fermion Hamiltonian we need to construct 
the space of antisymmetric mass matrix M that satisfy 

M 2 = -1. (78) 

The mass matrices A always anti commute with 7- 
matrices 7,, i = 1,2, ...,d. When the mass matrices M 
has some symmeties, then the mass matrices satisfy more 
linear conditions. Let us assume that all those conditions 
can be expressed in the following form 

Mpi = -piM, pipj = -pjpi, 
MUj = UjM, U lPl = ptUi, (79) 

where pi and XJj are real matrices labeled by i and J, 
and 71, ...,7d are included in p^s. If we have another 
symmetry condition W such that MW = —WM and 
Wpi = Pi W for a partcular i , Then U = Wpi will 
commute with M and pi, and will be part of Uj. 

Uj will form some algebra. Let us use a to label the 
irriducible representations of the algebra. Then the one 
fermion Hilbcrt space has a form H = ©a^a^"^, where 
the space Hq forms the a th irreducible representations 
of the algebra. For such a decomposition of the Hilbert 
space, M has the following block diagonal form 

M = ® a (M a ®I a ), (80) 

where I a acts within W® a as an identity operator, and 
M a acts within H. a . piS have a similar form 

p t = ® a (p? ®I a ), (81) 

where pf act within H a . So within the Hilbert space T~L ai 
we have 

Pi = ~Pl M ) Pi P 3 = ~Pj Pi ■ i^ 2 ) 

What we are trying to do in this paper is actually to con- 
tract the space of M a matrices that satisfy the condition 
eqn. (82). 

If fermions only form one irreducible representation of 
the Ui algebra, then the classifying space of M a and M 
will be the same. The results of this paper (such as Tables 
II, III) are obtained under such an assumption. 

If fermions form n distinct irreducible representations 
of the Uj algebra, then the classifying space of M will be 
R n , where R is the classifying of M a constructed in this 
paper. Note that the classifying spaces of M a are the 
same for different irreducible representations and hence 
R is independent of a. So if M a 's are classified by Z^, 
k = 1,2, or 00, then M's are classified by Z£. 

To illustrate the above result, let us use the symme- 
try G+(C) = Z2 x Z( as an example. If the fermions 
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form one irreducible representation of , for exam- 
ple, CciC'" 1 = — Cj, then the non-interacting symmetric 
gapped phases are classified by 



d : 

gapped phases 



1 2 3 4 5 6 7 
Z 2 Z 2 Z Z 



(83) 



which is the result in Table III. If the fermions form both 
the irreducible representations of Z§ , Cci+C^ 1 = +q + 



and Ccj-C 



(ie one type of fermions carries 



Z^-charge and another type of fermions carries Z 2 - 
charge 1), then the non-interacting symmetric gapped 
phases are classified by 



d; 12 3 4 5 

gapped phases : T\ T\ Z 2 



6 7 
Z 2 



(84) 



The four d = phases correspond to the ground state 
with even or odd Z 2 -charge-0 fermions and even or odd 
Z 2 -charge-l fermions. The four d — 1 phases correspond 
to the phases where the Z 2 -charge-0 fermions are in the 
trvial or non-trivial phases of Majorana chain and the Z 2 - 
charge-1 fermions are in the trvial or non-trivial phases 
of Majorana chain. The d = 2 phase labeled by two inte- 
gers (m, n) G Z 2 corresponds to the phase where the Z 2 - 
charge-0 fermions have m right moving Majorana chiral 
modes and the Z 2 -charge-l fermions have n right mov- 
ing Majorana chiral modes, (if m and/or n are negative, 
we then have the corresponding number of left moving 
Majorana chiral modes.) 

Some of the above gapped phases have intrinsic 
fermionic topological orders. So only a subset of them 
are non-interacting fermionic SPT phases: 



d sp : 1 2 3 4 5 6 7 

SPT phases : Z 2 Z 2 Z Z 



(85) 



The two d S p = phases correspond to the ground states 
with even numbers of fermions and or 1 Z 2 -charges. 
The two d sp = 1 phases correspond to the phases where 
the Z 2 -charge-0 fermions and the Z 2 -charge-l fermions 
are both in the trvial or non-trivial phases of Majorana 
chain. The d sp = 2 phase labeled by one integers n G Z 
corresponds to the phase where the Z 2 -charge-0 fermions 
have n right moving Majorana chiral modes and the Z 2 - 
charge-1 fermions have n left moving Majorana chiral 
modes. 



IV. CLASSIFICATION WITH TRANSLATION 
SYMMETRY 

For the free fermion systems with certain internal sym- 
metry Gf, we have shown that their gapped phases are 
classified by ■na{R PG -d) or 7To(C PG _d) in d dimensions, 
where the value of pa is determined from the full sym- 
metry group Gf. We know that ttq is an Abelian group 
with commuting group multiplication "+" : a, b G 7r im- 
plies that a + b G ttq. The "+" operation has a physical 



meaning. If two "bulk" gapped fermion systems are la- 
beled by a and b in ttq, then stacking the two systems 
on together will give us a new "bulk" gapped fermion 
system labeled by a + b G ttq- 

We note that the classification by 7ro(i? PG -d) or 
TTo(Cp G —d) is obtained by assume there is no translation 
symmetry. In the presence of translation symmetry, the 
gapped phases are classified differently. 32-34,42 However, 
the new classification can be obtained from 7ro(i? PG _(i). 
For the free fermion systems with internal symmetry Gf 
and translation symmetry, their gapped phases are clas- 
sified by 25 



(86) 



fe=0 



where (f) is the binomial coefficient. The above is for the 
real classes. For the complex classes, we have a similar 
classification: 



f[[Mc PG - d+k )p. 



(87) 



fe=0 



Such a result is obtained by stacking the lower dimen- 
sional topological phases to obtain higher dimensional 
ones. For 1-dimensional free fermion systems with inter- 
nal symmetry Gf, their gapped phases are classified by 
7r (i? PG _i). We can also have a 0-dimensional gapped 
phase on each unit cell of the 1-dimensional system 
if there is a translation symmetry. The 0-dimcnsional 
gapped phases are classified by itq(R pg ). Thus the com- 
bined gapped phases (with translation symmetry) are 
classified by 7ro(-R Po -i) x 7ro(-Rp G )- In 2-dimensions, the 
gapped phases are classified by 7To(i? PG - 2 ). The gapped 
phases on each unit cell are classified by ttq(R pg ). Now 
we can also have 1-dimensional gapped phases on the 
lines in x-direction, which are classified by 7To(i? PG _i). 
The same thing for the lines in y-direction. So the 
combined gapped phases (with translation symmetry) 
are classified by 7r (i? PG _ 2 ) x [7r (-R PG _i)] 2 x 7r (i? PG ). 
In three dimensions, the translation symmetric gapped 
phases are classified by 7r (i? PG _ 3 ) x [7r (_R PG _ 2 )] 3 x 



[n Q (R. 



X TT {R pG ). 



V. DEFECTS IN d-DIMENSIONAL GAPPED 
FREE FERMION PHASES WITH SYMMETRY G f 

For the cZ-dimensional free fermion systems with in- 
ternal symmetry G/, we have shown that their "bulk" 
gapped Hamiltonians (or the mass matrices) form a space 
R PG -d or C PG _rf. (More precisely, the gapped Hamiltoni- 
ans or the mass matrices form a space that is homotopi- 
cally equivalent to R PG -d or C PG -d)- From the space 
R PG -d or C PG ~d, we find that the point defects that 
have the symmetry Gf are classified by ^d-i{R PG -d) or 
7i"d-i(C PG _d). 
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Symmetry 




pa 


point 
defect 


line 
defect 


example phases 


U(l) 
G-(C) 


U(l + m) -, 

U(l)xU(m) A L 








Z 


(Chern) supercond. 

with collinear 
insulator gpin order 


Gf_(U,T) 
G±_(T, C) 

gX_(t, C) 


U(n) 


1 


z 





supercond. w/ real pairing 
and S z conserving 
spin-orbital coupling 



TABLE VI: (Color online) Classification of point defects and line defects that have some symmetries in gapped phases of 
non-interacting fermions. "0" means that there is no non-trivial topological defects. Z means that topological non-trivial 
defects plus the topological trivial defect are labeled by the elements in Z. Non-trivial topological defects have protected 
gapless excitations, while trivial topological defects have no protected gapless excitations. 



Symm. 


G~(U,T) 
G-_{T,C) 


G+(T) 
G+ + (T, C) 
G~-{U,T, C) 
G- + (U,T,C) 
G±-(U, T, C) 
G++(U,T, C) 


"none" 
G+(07) 
G- + (T,C) 
GZ+(T, C) 

G+(C) 
G + {U,C) 


G-(T) 
G± + (T,C) 
GZ+(U,T, C) 
GZ\(U, T,C) 
G\Z{U, T,C) 
G±+(U,T, C) 


GZ(U, T) 
GI_(T, C) 


GZZ(U, T,C) 
GZt(U, T, C) 
G±Z(U, T,C) 
G+±(U,T, C) 


G-(U,C) 
SU(2) 


GZZ(U, T, C) 
GZt(U,T, C) 
G%Z{U,T, C) 
G±t(U, T, C) 
G[SU(2),T] 


R PG 


0(l+m) T 


0{n) 


0(2n) 


U(2n) 


Sp(l+m) 


Sp{n) 


Sp(n) 


U(n) 


0(l)xO(m) A *- 


U(n) 


Sp(n) 


Sp(l)xSp(m) A "~ 


C7(n) 


PG- 





1 


2 


3 


4 


5 


6 


7 


point 
defect 





Z 


z 2 


z 2 





Z 








line 
defect 








Z 


z 2 


z 2 





Z 





example 
phases 


insulator 
w/ coplanar 
spin order 


supercond. 
w/ coplanar 
spin order 


supercond. 


supercond. 
w/ time 
reversal 


insulator 
w/ time 
reversal 


insulator 
w/ time 
reversal and 
inter-sublattice 
hopping 


spin 
singlet 
supercond. 


spin 
singlet 
supercond. 
w/ time 
reversal 



TABLE VII: (Color online) Classification of point defects and line defects that have some symmetries in gapped phases of non- 
interacting fermions. "0" means that there is no non-trivial topological defects. Z n means that topological non-trivial defects 
plus the topological trivial defect are labeled by the elements in Z„. Non-trivial defects have protected gapless excitations in 
them. 



Physically, there is another way to classify point de- 
fects: we can simply add a segment of ID "bulk" gapped 
free fcrmion hopping system with the same symmetry to 
the d-dimensional system. Since the translation symme- 
try is not required, the new c£-dimensional system still 
belong to the same symmetry class. There are finite 
bulk gap away from the two ends of the added ID seg- 
ment. So the new d-dimensional system may contain two 
non-trivial defects. The classes of defect is classified by 
the classes of the added ID "bulk" gapped free fermion 
hopping system. So we find that the point defects that 
have the symmetry Gf are also classified by Tro{Rp G -i) 
or 7r (C PG _i). 

Similarly, the line defects that have the symmetry Gf 
are classified by TT d ^ 2 {R PG -d) or r Kd—2(Cp G —d)- Again, 
we can also create line defects by adding a disk of 2D 
"bulk" gapped free fermion hopping system to the origi- 
nal ci-dimensional system. This way, we find that the line 
defects that have the symmetry Gf are also classified by 



7r (i?p G _ 2 ) or n (C PG - 2 )- 

In general, the defects with dimension do are classified 
by nd-d -i(Rp G -d) or Kd-do-i{ C PG-d), or equivalently 
by iro(R PG -d -i) or Tro(C PG -d -i)- In order for the above 
physical picture to be consistent, we require that 

ir n (R PG ) = n (R PG+n ), ir n (C PG ) = Tr (Cp G+n ). (88) 

The classifying spaces indeed satisfy the above highly 
non-trivial relation. This is the Bott periodicity theorem. 
The theorem is obtained by the following observation: 
the space C p +i can be viewed (in a homotopic sense) as 
nCp - the space of loops in C p . So we have ni(C p ) — 
7To(fiC p ) = 7To(C p +i). Similarly, the space R p +i can be 
viewed (in a homotopic sense) as £lR p - the space of loops 
in Rp. So we have iri(R p ) = ttq(URp) = ttq(R p+ i). As 
a result of Bott periodicity theorem, the classification of 
defects is independent of spatial dimensions. It only de- 
pends on the dimension and the symmetry of the defects. 
If the defects lower the symmetry, then, we should use 
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the reduced symmetry to classify the defects. In Tables 
VI and VII, we list the classifications of those symmetric 
point and line defects for gapped free fermion systems 
with various symmetries. We would like to point out 
that the line defects classified by Z in superconductors 
without symmetry do not correspond to the vortex lines 
(which usually belong to the trivial class under our clas- 
sification). The non-trivial line defect here should carry 
chiral modes that only move in one direction along the 
defect line. In general, non-trivial defects have protected 
gapless excitations in them. 



reversal T, and/or charge conjugation C symmetry. We 
show that each symmetry group Gf is associated with a 
classifying space C PG or R PG (see Tables II and III). We 
classify d-dimensional gapped phases of free fermion sys- 
tems that have those full symmetry groups. We find that 
the different gapped phases are described by 7ro(C PG _,j) 
or Tra(Rp G -d)- Those results, obtained using the K- 
theory approach, generalize the results in Ref. 25,28 



VI. SUMMERY 

In this paper, we study different possible full symme- 
try groups Gf of fermion systems that contain U(l), time 
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